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1 Introduction 



The aim of this review is to study gluon scattering amplitudes of four dimensional 
planar maximally super-symmetric Yang-Mills (MSYM). We hope that the study of 
such amplitudes would teach us something about scattering amplitudes of QCD, but 
at the same time they are much more tractable. The reason for such tractability is 
twofold. On one hand, perturbative computations are much simpler than in QCD, 
due to the high degree of symmetry. In fact enormous progress has been made in 
the last few years. On the other hand, the strong coupling regime of the theory can 
be studied by means of the AdS/CFT duality, by studying a weakly coupled string 
sigma-model. 

In this review we focus on how to use the AdS/CFT duality in order to compute 
gluon scattering amplitudes of planar MSYM at strong coupling, referring the reader 
to I V.l , V.2| for details on the perturbative side of the computation . In section two 
we set up the problem of computing scattering amplitudes at strong coupling. The 
problem boils down to the computation of the area of certain minimal surfaces in AdS. 
For the particular case of four gluons, such surface, and its area, can be explicitly 
computed. Furthermore, the strong coupling computation hints at some symmetries 
that actually appear to be symmetries at all values of the coupling. This is briefly 
reviewed at the end of section two. In section three we focus on the mathematical 
problem of computing the area of minimal surfaces in AdS. The integrability of the 
model allows the introduction of a spectral parameter. By studying the problem as a 
function of the spectral parameter we are able to give a solution in the form of a set of 
integral equations. These equations have the precise form of thermodynamic Bethe 
ansatz (TBA) equations. The area turns out to coincide with the free energy of such 
TBA system. Finally, In section four, we end up with some conclusions and a list of 
open problems. 

2 Gluon scattering amplitudes at strong coupling 

Four dimensional MSYM, the theory whose amplitudes we want to consider, turns 
out to be dual to type IIB string theory on AdS 5 x S^. This duality receives the name 
of AdS/CFT duality [ 1 1 and is the main focus of this review. A remarkable feature 
of this duality is that it allows to compute certain observables of MSYM at strong 
coupling by doing geometrical computations on AdS. A well known example is the 
computation of the expectation value of super-symmetric Wilson loops, which re- 
duces to a minimal area problem [2, 3 1. In this section we will show that this is also 
the case for the computation of scattering amplitudes at strong coupling ! ^ 

As in the gauge theory, we will need to introduce a regulator in order to define 
properly scattering amplitudes. In order to set-up our computation we introduce a 

'in this section, we follow closely |4|, to which we refer the reader for the details. 
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D-brane as IR regulator, as we explain in detail below. Another convenient regula- 
tor is the strong coupling/super-gravity analog of dimensional regularization. This 
regulator will be used in order to compare our results with expectations from the 
perturbative side. 

2.1 Set-up of the computation 

In order to set up the computation at strong coupling, it is convenient to introduce 
a regularization as follows. We start from a.U{N + k) theory, with k <^N, and then 
consider a vacuum breaking the symmetry to U (N) x U (k) by giving to a scalar 
field a vacuum expectation value mjR which plays the role of an infrared cut-off 
When we take the 't Hooft limit we keep k fixed, so that the low energy U (k) theory 
becomes free. We then scatter gluons of this U {k) theory. We are interested in the 
regime where all kinematic invariants are much larger than the IR cut-off, Sij » mjjf. 
It turns out that the leading exponential behavior at strong coupling can be captured 
simply by considering k = I. At strong coupling this corresponds to consider a D3- 
brane localized in the radial direction. More precisely, we start with the AdS^ metric 
written in Poincare coordinates 

^s^=R2 d4+i+dz^ (2.1) 

and place a D3-brane at some fixed large value of z = zir and extending along the 
;c3+i coordinates. The asymptotic states are open strings that end on that D-brane. 
We then consider the scattering of these open strings, that will have the interpretation 
of the gluons that we are scattering. 

The proper momentum of the strings is kpr = kzm/R, where k is the momentum 
conjugate to X3+1, plays the role of gauge theory momentum and will be kept fixed 
as we take away the IR cut-off, zm — )• °°- Therefore, due to the warping of the metric, 
the proper momentum is very large, so we are considering the scattering of strings at 
fixed angle and with very large momentum. 

Amplitudes in such regime were studied in flat space by Gross and Mende |6|. 
The key feature of their computation is that the amplitude is dominated by a saddle 
point of the classical action. In our case we need to consider classical strings on 
AdS. Hence, we need to consider a world-sheet with the topology of a disk with 
vertex operator insertions on its boundary, which correspond to the external states 
(see fig. 1). A disk amplitude with a fixed ordering of the open string vertex operators 
corresponds to a given color ordered amplitude. 

What are the boundary conditions for such world-sheet? since the open strings 
are attached to the D-brane, z = Zm at the boundary. Furthermore, in the vicinity 

^See |5 1 for perturbative computations using this regulator. 
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Figure 1: World-sheet corresponding to the scattering of four open strings. 

In the figure on the left we see four open strings ending on the IR D-brane, 
the world- sheet has then the topology of a disk, shown on the right, with four 
vertex operator insertions. 

of a vertex operator, the momentum of the external state should fix the form of the 
solution. 

In order to state more simply the boundary conditions for the world-sheet, it 
is convenient to describe the solution in terms of T-dual coordinates y^, defined as 
follows 



The presence of the / is due to the fact that we are considering a Euclidean world-sheet 
in Minkowski space-time. Note that we do not T-dualize along the radial direction. 
After defining r = R^/z the dual metric takes the form 



Note that this metric is equivalent to the same AdSs metric we started with! A crucial 
difference is that now, in terms of the dual coordinates, the boundary of the world- 
sheet is located at r = /zir, which is very small. Furthermore, the T— duality 
we performed interchanges Neumann by Dirichlet boundary conditions. This means 
that the boundary of the world-sheet sits at a fixed point in the space of the dual 
coordinates. When a vertex operator with momentum is inserted, the location of 
such point gets shifted by an amount proportional to = Ink^. 

Summarizing, the boundary of the world-sheet is located at r = /zm and is a 
particular line constructed as follows 

• For each particle of momentum k^, draw a segment joining two points sepa- 
rated by Aj'^ =2nk^. 



ds^ = {z)dx^dx^ + ... day^ = iw^ {z)£ap dpx' 



(2.2) 



ds^ = R 



^2dy^dy^+dr^ 
r2 



(2.3) 
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• Concatenate the segments according to the insertions on the disk. 

Since gluons are massless, the segments are light-like. Furthermore, due to mo- 
mentum conservation, the segments form a closed polygon. The world-sheet, when 
expressed in T-dual coordinates, will then end in such sequence of light-like segments 
(see fig. 2) located at r = R^/zm- 




Figure 2: Comparison of the world sheet in original and T-dual coordinates. 
The hyperplane on the picture to the right should not be interpreted as a D- 
brane, but rather as a radial slice where the boundary of the world-sheet it 
located. 

As we take away the IR cut-off, z.ir — )■ °°, the boundary of the world-sheet moves 
towards the boundary of the T-dual metric, at r = 0. This computation would then 
be formally equivalent to the computation of the expectation value of a Wilson loop 
given by a sequence of light- like segments at strong coupling [2,3 |. ^ 

Our prescription is that the leading exponential behavior of the A'^— point scatter- 
ing amplitude is given by the area A of the minimal surface that ends on a sequence 
of light-like segments on the boundary 

e-^^^''-'^^ (2.4) 

An important comment is in order. Note that the strong coupling computation is blind 
to the type or polarization of the external particles. Such information will contribute 
to prefactors in (2.4) and will be subleading in a expansion, relative to the 

leading exponential term^ These differences should be visible once we consider 
quantum corrections to the classical area. This is still an open problem. 

^As explained in detail in | V.2|, this remarkable duality between Wilson loops and scattering ampli- 
tudes was also observed in perturbative computations. 

'^As discussed in more detail in |7|, when computing the disk amplitude in the saddle point approx- 
imation, one can neglect the polarization of the gluon vertex operators. 
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Some generalizations to the above picture were developed. In [ 8 ,9 1 finite temper- 
ature was introduced while in |10| the authors considered solutions with non trivial 
motion on the S^. Finally, in 1 1 1-13], the scattering of quarks at strong coupling was 
considered. Unfortunately, due to space limitations, we wont discuss these interesting 
developments here, but refer the reader to the original literature. 

We have then reduced the problem of computing scattering amplitudes at strong 
coupling to the problem of finding minimal surfaces in AdS. In the following we will 
show that such surface can be found for the particular case of the scattering of four 
gluons. To find and understand this solution in detail will be quite instructive. Then, 
in the next section, we will use the integrability of the problem in order to give a 
general solution, for any number of gluons, in the form of a set of integral equations. 

2.2 Scattering of four gluons 

Consider the scattering of two particles into two particles, + ^3 — )■ ^2 + ^4 and 
define the usual Mandelstam variables 

S = -{ki+k2)\ t = -{k2+hf (2.5) 

According to our prescription we need to find the minimal surface ending in the 
following light-like polygon 



Figure 3: Polygon corresponding to the scattering of four gluons 

In order to write the Nambu-Goto action it is convenient to use Poincare coordi- 
nates {r,yQ,y\,y2), setting = and parametrize the surface by its projection to the 
{y\,y2) plane. In this case we obtain an action for two fields, r and t, living in the 
space parametrized by y\ and y2 

. f , , yi + [djrf - {djyof - {d,rd2yo - d2rd,yoY 

Ik I -^^^^ ^ ' 

where r = r{yi,y2) and = dy.r, etc. Our aim is to find a solution to the classical 
equations of motion with the appropriate boundary conditions. Let us consider first 
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the case s = t, where the projection of the polygon Unes to the {y\,y2) plane is a 
square. By scale invariance we can choose the edges of the square to lie at y\,y2 = 
zb 1 . The boundary conditions are then given by 



r(±l,3;2)=r(3;i,±l)=0, 3'o(±l,3'2) = ±3^2, 3'o(3'i , ±1) = ±^1 (2-7) 

In [14 1 the solution corresponding to a single cusp was considered. One can make 
educated guesses using such solution as a guidance and propose 

3'o(3'i,3'2) =3'i3'2, r{yi,y2) = ^^(1 -yfj{'^ -yl) (2.8) 

Remarkably this turns out to be a solution of the equations of motion! This is the 
solution for the case s = t, how can we obtain the most general solution? 

The ia&lAdSs space has a 50(2,4) group of isometrics. This symmetry is some- 
times referred to as "dual conformal symmetry" and should not be confused with the 
original 50(2,4) symmetry associated to the original AJS' space. This dual symme- 
try can be used in order to map the particular solution we have just found to the most 
general solution with four edges, in particular with s^t. The general solution can be 
conveniently written as 

n a\J l+b^ sinhwi sinhM2 

3^0 = — r r , , ■ . r-r — (2.9) 



cosh M 1 cosh U2+b sinh u \ sinh U2 cosh m i cosh U2 + b sinh u \ sinh U2 

asinhwi coshM2 acoshMi sinhM2 f2 10) 

cosh Ml cosh M2 +^ sinh Ml sinh M2 ' "^^ cosh m i cosh M2 + ^ sinh mi sinh M2 

where mi 2 parametrize the world-sheet and we have written the surface as a solu- 
tion to the equations of motion in conformal gauge, a and b encode the kinematical 
information of the scattering as follows 



According to the prescription, we should now plug the classical solution into the 
classical action to compute the area and obtain the four point scattering amplitude 
at strong coupling. However, in doing so, we obtain a divergent answer. That is of 
course the case, since we have taken the IR regulator away. In order to obtain a finite 
answer we need to reintroduce a regulator. Since we want to compare our results to 
field theory expectations, it is convenient to introduce the strong coupling analog of 
dimensional regularization. 

Gauge theory amplitudes are regularized by considering the theory in D = 4 — 2£ 
dimensions. More precisely, one starts with ^ = 1 in ten dimensions and then di- 
mensionally reduce to 4 — 2£ dimensions. For integer 2e this is precisely the low 



7 



energy theory living on a D/j— brane, where /? = 3 — 2e. We regularize the ampli- 
tudes at strong coupling by considering the gravity dual of these theories and then 
analytically continuing in e. The string frame metric is 

Following the steps described above, we are led to the following action 

Where ^e=o is the Lagrangian density for AdSs. The presence of the factor 
will have two important effects. On one hand, previously divergent integrals will now 
converge (if e < 0). On the other hand, the equations of motion will now depend on 
£ and we were not able to compute the full solution for arbitrary e. However, we are 
interested in computing the amplitude only up to finite terms as we take e — )• 0. In 
that case, it turns out to be sufficient to plug the original solution into the £-deformed 
action ^ After performing the integrals and expanding in powers of e we get the final 
answer 



■e-^'', ,-— A = iSjn, + — [log-] +C (2.14) 




- 2Sdiy^s + 2Sdiv,t 



2e 



This answer has the correct general structure (see e.g. [V. 1 , V.2 1) from field theory 
expectations. Furthermore, once we use the strong coupling behavior for the cusp 
anomalous dimension |15|, /(A) = ^ + ... we see that the leading divergence, as 
well as the finite piece, have not only the correct kinematic al dependence but also the 
correct overall coefficient in order to match the BDS ansatz |V.2, 16 1. As we will see 
shortly, this result is actually a consequence of the symmetries of the problem! 



2.3 r— duality and dual conformal symmetry at strong coupling 

An important ingredient of the previous computation was the existence of a dual 
50(2,4) symmetry ^, associated to the isometry group of the dual AdSs space. This 

^Up to a contribution from the regions close to the cusps that adds an unimportant additional constant 
term. 

''Actually, this symmetry was first noticed in perturbative computations 1 17 1 and then independently 
in the strong coupling computation described here. 
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symmetry allowed the construction of new solutions and fixed somehow the finite 
piece of the scattering amplitude. ~' 

To a symmetry we associate a Ward identity and in particular dual conformal 
symmetry will impose some constraints on the amplitudes. Quite remarkably, this 
duality was also (actually before!) observed at week coupling and is by now believed 
to be a duality of scattering amplitudes at all values of the coupling. You can see [ V.2 | 
for a detailed account of this symmetry and the constraints it imposes on the ampli- 
tudes.^ Here we will just mention that dual conformal symmetry fixes the answer 
for the four-point function to have the form (2.14 ), actually, to all values of the cou- 
pling! and hence its agreement with the BDS ansatz. Furthermore, dual conformal 
symmetry does not fix the answer for the scattering of more than six gluons, hence, 
in general, the answer deviates from the BDS ansatz. The need for such a deviation, 
usually called remainder function, was established in [20,21]. See |V.2| for more 
details. 

In the last section we have seen that existence of a dual AdS space, is related to 
the fact that AdSs goes to itself after a sequence of four T-dualities, followed by the 
inversion of the radial coordinate, see (2.1 ) vs (2.3). This set of T-dualities, however, 
does not leave the full AdS=, x sigma model invariant. For instance, Buscher rules 
for T-dualities 1 22] imply a shift on the dilaton of the form 

0^0+4xlogz (2.15) 

where z is the radial coordinate of the original metric (2.1 ). The factor of 4 is due 
to the fact that we are making four 7— dualities. In addition to the usual, "bosonic", 
T— dualities, one can introduce a fermionic T— duality |23,24|. This duality is a non 
local redefinition of the fermionic world-sheet fields, very much like the bosonic T- 
duality is a redefinition of the bosonic fields. These T-dualities change the fields of 
the sigma model according to precise rules. For instance, each fermionic T— duahty 
shifts the dilaton by an amount 

0^0-^xlogz (2.16) 

We see that by doing eight fermionic T— dualities we can undo the shift (2. 15 ) on the 
dilaton. Actually, one can check that a combination of the four bosonic T— dualities 
plus eight fermionic T— dualities maps the full sigma model to itself! Note also that 
this argument does not depend on the value of the coupling. One of the implications is 

^Naively, this conformal symmetry would imply that the amplitude is independent of s and t, since 
they can be sent to arbitrary values by a dual conformal transformation. The whole dependence on s 
and t arises due to the necessity of introducing an IR regulator. However after keeping track of the 
dependence on the IR regulator, the amplitude is still determined by the dual conformal symmetry. 
Hence, this regulator breaks the dual conformal symmetry, but in a controlled way! 

^These constraints have also been derived at strong coupling 1 18 19|. 
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that the dual model has the same conformal symmetry group as the original, helping 
to understand the origin of dual conformal symmetry. Actually, as the construction 
suggests, dual conformal symmetry extends to a full dual super conformal symmetry. 
In addition, one has a map between the full set of conserved charges of the two 
models, in such a way that some of the local charge of one model are mapped to non 
local charges of the dual model, and viceversa, see for instance |23,24|. 

The structure of dual super conformal symmetry was also seen at weak coupling 
and is explained in detail in | V.2 |, for which we refer the reader for more details. 

3 Minimal surfaces on AdS 

In the previous section we have seem how the problem of computing gluon scatter- 
ing amplitudes at strong coupling reduces to the computation of the area of certain 
minimal surfaces in AdS. In this section we show how the integrability of the system 
can be used in order to give a solution to the problem, in the form of a set of inte- 
gral equations. We will follow closely |25-28|, see also |29,30|, to which we refer 
the readers for the details. ^ For this review, we will focus mostly on a particular 
kinematic configuration, in which the minimal surfaces are actually embedded into 
an AdS^ subspace of the full AdSs. However, the full problem has been solved and it 
will be briefly mentioned at the end of the section. 

The mathematical problem is to find the area of the minimal surface ending on the 
boundary of AdS at a given polygon of light-like edges. The polygon is parametrized 
by the location of its cusps Xi, which are null separated, namely xf^^^ = 0. 

We will focus on certain regularized area that is invariant under conformal trans- 

formations. As such, it will depend only on cross-ratios, of the form Xuki = "i-r- 

Given the cross-ratios, we want to compute the area as a function of those. The full 
problem involves minimal surfaces on AdS 5, in which case there are 3N — 15 in- 
dependent cross-ratios, where N is the number of cusps/gluons. We will restrict to 
special kinematical configurations in which the minimal surfaces involved are em- 
bedded in AdS3. In this case, we have N — 6 independent cross-ratios and the 
polygon is a zig-zaged polygon living in one plus one dimensions, which correspond 
to the boundary of AdSs, see figure 4 

Since we want a closed contour, and we are in 1 + 1 dimensions, we can consider 
only polygons with an even number of sides, hence N = 2n. As one can see in figure 

^Some of the key ideas used below may be found in relation to the study of wall crossing |31 ,32|. 
Actually, the method of the first paper in 1 3 1 , 32 1 where instrumental in deriving the expression for the 
eight gluon amplitude at strong coupling in 1 25 1. 

"'For the general scattering in four dimensions we have 4A' coordinates, minus A', since the distance 
between consecutive points has to be light-like, minus 15, that is the dimension of the conformal group 
50(2,4). In the case of AdSj, we have 2N — N minus 6, which is the dimension of 50(2, 2). 
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Figure 4: A zig-zaged null polygon in 1 + 1 dimensions is parametrized by n 
xf' coordinates and n xj coordinates. If you want a closed polygon, you can 
fold the figure in a cylinder 



4 the contour is parametrized by n coordinates x[ and n coordinates . With each 

set of coordinates we can form « — 3 invariant cross-ratios, of the form y. , , = -^-1. 

In order to consider minimal surfaces in AdS^ we need to consider the world-sheet 
of classical strings on A<i53. This is the subject of the following subsection. 



3.1 Strings on AdSs 

Classical strings on AdS^, can be described in terms of embedding coordinates, where 
AdS3 is the following surface embedded in 7?^'^ 

Y.Y = -Y^^-Y^ + Y^ + Y^ = -\ (3.1) 

we take the world-sheet to be the whole complex plane. Since we are interested in 
classical solutions, the fields have to satisfy the conformal gauge equations of motion 
and the Virasoro constraints 

ddY -{dY.dY)Y = 0, dY.dY = dY.dY = (3.2) 

where dY = d^Y, etc. An efficient way to focus only in the physical degrees of 
freedom, similar to fixing light-cone gauge, is by performing the so-called Pohlmeyer 
kind of reduction, see for instance |33 -35|, and consider the "reduced" fields 

a = logdY.dY, p^ = d^Y.d^Y (3.3) 

As a consequence of the equations of motion and Virasoro constraints, p can be 
seen to be a holomorphic function, p = p{z) while a{z,z) can be seen to satisfy a 
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generalized version of the Sinh-Gordon equation 



dda-e''+p{z)p{z)e-'' =0 



(3.4) 



From the definition of the reduced fields, it is clear that they are invariant under space- 
time conformal transformation. This means that they describe only the essential part 
of the problem, without redundancies. 

Before proceeding, let us make the following remark. Since p{z) is a holomorphic 
function, it is possible to make a change of coordinates from the z — plane to the 
w — plane, where dw = a/ p{z)dz. In the w— plane, after a simple field redefinition, 
the generalized Sinh-Gordon equation takes the usual form 



It would seem that we got rid of all the information on p{z). However, this is not the 
case, since the w— plane will have in general a complicated structure (for instance, 
it will have a branch cuts, etc, depending on p{z)). So, we can choose between a 
complicated equation on the complex plane, or a simple equation on a more compli- 
cated space. Depending which questions we want to answer, one description may be 
more convenient than another. Finally, we are interested in the area of the classical 
World- sheet. Written in terms of the reduced fields it becomes 



3.2 Classical solutions corresponding to minimal surfaces ending on 
null polygons 

What are the properties of the holomorphic function p{z) and a{z,z) for solutions 
corresponding to minimal surfaces ending on null polygons? In order to answer this 
question we can start by considering the four cusps solution found in the previous 
section and perform the Pohlmeyer reduction. We find 



Hence, the four cusps solution simply correspond to the vacuum solution of the Sinh- 
Gordon equation! What about solutions with a higher number of cusps? First of all 
we propose that the field a is regular everywhere, since we are looking for smooth 
space-like solutions. Second, we expect a general solution to be similar to the four 
cusps solution when approaching the boundary, so we expect that a — > as |z| be- 
comes large. 



a = a + - log pp — )• dwdwOl = 2 sinh a 



(3.5) 




(3.6) 



;7(z) = l, a = a = 



(3.7) 
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Finally, if we are interested on a minimal surface ending on a polygon with 2n 
cusps, we propose p{z) to be a polynomial of degree « — 

p{z)=z:'-^ + Cn-Az"-^ + ...+ciz + co (3.8) 

we have used rescalings and translations in order to set the coefficients of z"^^ and 
z"^^ to one and zero respectively. Such polynomial contains « — 3 complex coeffi- 
cients, or 2?i — 6 real coefficients, which exactly agrees with the amount of expected 
independent cross ratios for a polygon with 2n cusps ! 

Summarizing: minimal surfaces ending on a light-like polygon with 2n cusps 
correspond to a holomorphic polynomial of degree n — 2 and a field a satisfying the 
Sinh-Gordon equations and with boundary conditions such that it decays at infinity 
and diverges logarithmically at the zeroes of p{z), which amounts to say that a is 
regular everywhere. 

Since a decays at infinity, the integral defining the area (3.6) diverges. We define 
a regularized area by subtracting the behavior at infinity 

Areg = /(^" - (3-9) 

As the reduced fields are invariant under space-time conformal transformations, the 
regularized area will be a function of the cross-ratios only. The computation of 
this regularized area is the main focus of the remaining of this review. 



3.2.1 Reconstructing the space-time solution and its behavior at infinity 

In the following we would to check that the world-sheet we are considering has the 
desired form. In particular, we would like to understand the shape, in space-time, 
of the boundary of our world-sheet. For that, we first review a general procedure to 
reconstruct the world-sheet from the reduced fields, and then study its boundary. 

Given an holomorphic function p{z) and a field a satisfying (3.4) it is possible to 
reconstruct a space-time solution satisfying (3.2), and (3.1). The procedure amounts 
to solve two auxiliary linear problems, which we denote as left and right 

{d+B'^)^if^,=0, {d + B'')Y! = (3.10) 

where the flat connections B^ ^^ are two by two matrices constructed from p{z) and 
a{z,z)- For instance 

' ' The motivation for this proposal, comes from the fact that an homogeneous polynomial of degree 
n — 2 possesses all the symmetries to correspond to a symmetric polygon of n edges. 

^^The full answer would include also the integral of the one we have subtracted. In order to compute 
it one would need to introduce a physical regulator and this part of the answer will not be conformal 
invariant. Anyway, its explicitly form can be worked out and turns out to be quite universal. In this 
review we will focus on the "interesting" part of the answer Areg. 
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we denote different components of the connections by B^^ and . On the other 
hand, the indices a and d in (3.10) denote independent solutions of the auxiliary 
linear problems. Each i//^ or is then a doublet. We denote the components of this 
doublet by Y^.a^ ^tc. 

Given the solutions of these two auxiliary linear problems, one can show that the 
space-time coordinates are simply given by 

One can show that Y constructed this way satisfies all the required properties. If we 
see 1//^ and two by two matrices, then the space-time coordinates would be 

given by Y = {y^Yy^. On the other hand, note that given a solution to the left 
problem, i//^, then Xj/^U^ is an equally good solution, and the same happens with 
the right problem. Hence, given Y, we obtain a family of space-time solutions UYV. 
These are nothing but the space-time conformal transformations. 

Now we would like to understand the behavior of the solutions of the linear aux- 
iliary problems for very large values of |z|, or |w|. This will tell us the behavior of the 
world-sheet near its boundary. Let us start, by simplicity, with the case of a homoge- 
neous polynomial, p{z) = z"~^- Hence w ?a z"^^- As a result, as we go once around 
the z— plane, we go around the w— plane n/2 times. 

Due to the boundary conditions for the reduced fields, the flat connections B^'^ 
drastically simplify at infinity and we can solve the auxiliary linear problems. A 
general solution will be of the form 

W! ^ 4 Q +dT (^^ e-'-^ (3.13) 

The w-plane is naturally divided into quadrants, see figure 5 In each quadrant one 
of the two solutions of each problem (left and right) dominates. For instance, in the 
upper right quadrant, the solution proportional to c+ dominates in the left problem, 
while the solution proportional to dominates in the right problem. This means 
that for large values of \w\, the whole quadrant corresponds to a single point in the 
boundary, given by Yaa ~ (Large) x c+(i+. 

As we change quadrant, one and only one of the two dominant solutions change 
and we jump a light-like distance to the next cusp. In each quadrant/cusp we can 
write 

Ya,a-^ala (3.14) 

where A is given by the leading contribution to the left problem and A by the leading 
contribution to the right problem. As we change quadrant, one of the two solutions. 
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w-plane 




Figure 5: When looking at the left problem, (each sheet of) the w— plane is 
naturally divided into two parts, according to the sign of Re{w). In the same 
way, when looking at the right problem, the w— plane is naturally divided into 
two parts, according to the sign of Im(w). Hence, the w— plane is naturally 
divided into four quadrants. Large values of |w| in each of these angular sectors 
correspond to a cusp. 

A or A, changes. As we go around the w— plane n/2 times, we get the expected 2n 
cusps! 

In the general case in which the polynomial is not homogeneous, the picture is 
very much the same. In general, the degree of the polynomial determines the number 
of cusps, while the coefficients on the polynomial determine the shape of the polygon. 

Let us finish this section with a very important observation. Since we know 
that the classical equations describing strings on AdSs x are integrable, see for 
instance [37 1, we expect the present problem to be integrable as well. Indeed, it is 
possible to introduce a spectral parameter 

such that the flat connections are still flat, namely the satisfy dB^ — dB^ + [B^,B^ = 0, 
for all values of The introduction of the spectral parameter, also allows to study 
both linear problems in a unified manner. The left and right connections are just 
particular cases of the above flat connection, more precisely 

B(^ = 1)=B^, = (3.17) 
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The existence of the spectral parameter played a key role in the spectrum problem, see 
for instance [38 1. In that context, the key object are the eigenvalues of the monodromy 
matrix constructed out of the flat connection. In the present case the key object are 
certain cross-ratios constructed from the holonomy of the connections. 

3.3 Y-system for minimal surfaces 

Let us focus on the left problem. We see that each sheet on the w— plane is naturally 
divided into two sectors, one with Re{w) > and the other with Re{w) < 0. In each 
sector the small solution is well defined (up to a normalization constant). On the 
other hand, the large solution is not, as we can add to it a part of the small solution. 
Let us then introduce the following terminology: 

• The w— plane is divided into n sectors, since each sheet contains two sectors. 
We label this sectors by / = 0, — 1. 

• We call sf the small solution at the / — th sector. This is the solution with the 
fastest decay along the line in the center of the / — th sector, for increasing |w|. 

In order to understand why these small solutions are important, we need to introduce 
a new element. Given that our connections are SL{2) matrices, we can introduce a 
SL{2) invariant product 

Va^Vh^ £"'^ Wa^WpM = £ab (3. 18) 

The second equality corresponds to setting a normalization factor for our solutions. 
This can be done since one can check the above product is independent on the world- 
sheet coordinate z- 

As already seen, the location of the cusps is determined by the large solutions. 
The large component of a solution, on a given sector, can be extracted by using the 
small solution on such sector and the SL{2) invariant product just introduced, more 
precisely 

Wa^s\^Xi (3.19) 

where the symbol k, means up to factors that will cancel out in the final expression 
for the cross-ratios. How do we construct space-time cross-ratios? we have see that 
the location of the cusps is given by = A^A^. The space time cross-ratios involve 
distances like 

yi yj = £''beahY^.Yj.^ x'kj >< W >, < X'Xj >= e"^KXl (3.20) 
Given the normalization condition (3.18), one can easily show 

<X'V>^< Asf){\l/a As^) >=sfAs^ (3.21) 
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Which means that space-time cross-ratios can be constructed from inner products of 
the small solutions in the corresponding sectors ! more precisely 



X;,X 



Small solutions are defined up to a normalization constant. Note that such normal- 
ization constants cancel out when computing cross-ratios. 

The strategy we will follow is to introduce the spectral parameter as shown in 
the previous section and study the small solutions of the corresponding connection 

{d + B{Q)si{Q=0 (3.23) 

then, we can consider the cross-ratios as a function of such spectral parameter 

The physical cross-ratios are then obtained by setting the spectral parameter to ap- 
propriate values 

Zijkii^ = 1 ) = Xuki ' XijdC = i) = Xijki (3.25) 

A very important property of the flat connection B{Q is that it possesses a Z2 symme- 
try: B{e'^Q = GT,B{Qas, where 03 is the usual Pauli matrix. This symmetry allows 
to relate small solutions at different values of the spectral parameter, for instance 
^i+i{Q = <y3^i{^'^Q, and in particular, it implies 

SiAsj{e'''Q=Si+iAsj+,{Q (3.26) 

This identity is crucial in deriving the equations below. Besides, in order to simplify 
subsequent expressions, we will assume si A = 1. 

Now we have all the elements to derive the so called Hirota equations and the Y- 
system equations. The trick is to choose and as a complete basis of flat sections, 
and express two arbitrary consecutive small solutions and s^+i in terms of these 

Sk = {skAsi)so- {skAso)si (3.27) 
Sk+\ = {sk+iAsi)sQ-{sk+iAso)si (3.28) 

Next, use ( 3.26) in order to express every wedge as a wedge involving 50 and consider 
I = SkA Sk+i, we obtain 

- (sk-i Asoy^{sk+iAso) + {skAsoy+{skAso) = 1 (3.29) 
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where we have introduced the notation = /{e^'^^^Q, /++ = /{e'^^Q, etc. Let us 
introduce Tj^ = sqA s^+i (e^'^^+'^'^/^i^). In terms of these we obtain 

7;+r- = r,+i7:,_i + i (3.30) 

which has the form of the so called Hirota equations! from the definition of T^, we 
see that it is non trivial for s = 0,...,n — 2. The 7— system equations can be obtained 
by introducing = 7;_ir^+i 

Y+Y- = {\+Y,+i)i\+Y,-i) (3.31) 

Ys is non trivial for 5 = 1 — 3 . Note that this agrees with the amount of (complex) 
cross-ratios of our scattering problem. These are functional equations for Ys{Q and 
are valid for any value of ^. Note that they followed from a chain of rather trivial 
facts! 

One could reintroduce the normalized factors SjAsj^i and check that the F— functions 
are given by the usual cross-ratios introduced above. The physical cross-ratios, are 
then obtained by evaluating Ys{Q at = I and = /. 

Such equations are not the whole story. In particular, note that that the holomor- 
phic function p{z) does not enter at all in such equations! The point is the following. 
There are many solutions to such equations. The correct solution is then picked by 
specifying the analytic properties and boundary conditions of Ys{Q as we move on 
the ^— plane. This is how the information about the holomorphic polynomial enters 
and will be the subject of the following section. 

3.4 Integral equations 

In order to pick the appropriate solution to the 7— system equations (3.31) we need to 
specify the analytic properties of Fv(0- analyzing the auxiliary linear problems 
and the definition of Ys{Q one can show that Ys{Q are analytic away from = 0,°°. 
On the other hand. As — )• 0, °°, the fiat connection simplifies and the inverse problem 
can be solved by using a WKB approximation, where the role of h is played by ^ or 
l/(^ . By calling = e^, one can show that for large 6 the solution behaves as 

logF^ ^ -m.coshe + ... (3.32) 

where nta is given by the periods of pizY^^ along the cycles Ya, namely /Mq ^ — ^ p{z)dz. 
This is how the information of the polynomial piz) enters into the problem. This pe- 
riods are usually complex, and there are « — 3 of them, which exactly agrees with the 
quantity of expected cross-ratios. These should be seen as the boundary condi- 
tions for the above equations. 

'"'Note that even though we used the WKB approximation, this is the behavior of the exact solution. 
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The strategy now is well known from the study of integrable systems. We can 
combine the 7— system equations with the analytic properties and boundary condi- 
tions for the F— functions, in order to write a system of integral equations for them. 
The solutions to these integral equations will automatically satisfy the F— system 
equations and have the required boundary conditions. The system of integral equa- 
tions is given by 

I poo Y 

logF. = -»,.coshe + -y_^ ^^^^^^_^y og(l+F..,,(eO)(l+F.._,(eO) (3.33) 

The system of equations ( 3.33 ) has the form of TBA equations, that arise when study- 
ing integrable models in finite volume, see e.g. [36|. Even though, for the sake of 
clarity, some overall coefficients have been suppressed in the derivation of these 
equations, the final form of the equations is given with all the correct coefficients. 
From the TBA point of view, the parameters enter as masses. Once the masses 
are given, the solution of the above system is unique. The physical cross-ratios can 
the been read off by looking a.tYs{d) for appropriate values of 6. 

These integral equations can also be written in terms of the physical cross-ratios 
y+ = Fv(C = 1) ^iidj7 =Ys{^ = i) only, without resorting to the auxiliary parameters 
m.s. In order to achieve that, one simple evaluates (3.33 ) at the physical values of the 
spectral parameter in order to eliminate the masses, see |42|. We obtain 

logF, = cosh01og3;+-/sinheiogj; + J je'—^-^^^^^*^!-^ 

(3.34) 

Note that having solved (3.34), we could read off the masses from the asymptotic 
behavior of the solutions. 

How do we compute the regularized area, once we have solved the above sys- 
tem of integral equations? It turns out that the area can be written in terms of the 
F— functions in a very simple form 

/m 
J0^cosh01og(l+F,(0)) (3.35) 

In order to derive this expression, one expands the F functions a few orders for small 
and for large values of ^. The expansion coefficients are written in terms of period 
integrals that also appear in the expression for the area, see [28] for details. This 
expression has exactly the form of the free energy of the TBA system. , 

^'^Considering 4 = log(Kj/e^"'''^'"*''^), which is analytic in the strip Im{9) < Jz/2 and vanishes as 
9 approaches infinite. The integral equations can be obtained by convoluting the equation /+ + 1^ = 
log( 1 + 7,+ 1 ) ( 1 + y;_ 1 ) with the kernel in ( 3.33 ). 

'^It can be shown |42| that actually this area coincides with the extremum of the Yang- Yang func- 
tional for the modified TBA equations ( 3.34 ) 
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The strategy to solve the full problem is then clear. For a choice of the cross- 
ratios we solve the integral equations (3.34 ), and from their solution we compute the 
area (3.35). Hence, we have the area for these values of the cross-ratios. 

In this review we have treated in detail the case of minimal surfaces in AdSs. 
However, the general case of minimal surfaces in AdSs can also be solved [28]. Much 
of what we have said can be carried out for the general case. In this case we get a 
bigger system of F— functions, denoted by Ya.s, where a = 1,23 and 5=1, ...,N — 5. 
Note that their number equals the number of independent cross-ratios. Very much as 
before, one can obtain F— system equations, which supplemented with the appropri- 
ate boundary conditions can be written as a system of integral equations. Again, this 
system of equations has the form of a TBA system, and the regularized area coincides 
with the free energy of such system. 

These equations can be solved numerically, see for instance [28 1. On the other 
hand, it is very hard to find analytical solutions. However, some limits, for instance 
the so called small masses/CFT limit and the large masses limit , aie more tractable, 
see [39-41]. 

4 Conclusions 

We reviewed the computation of scattering amplitudes of planar maximally super- 
symmetric Yang-Mills at strong coupling. By using the AdS/CFT duality the prob- 
lem boils down to the computation of the area of certain minimal surfaces on AdS. 

Then we showed how the integrability of the model can be efficiently used in 
order to give an answer for the problem in terms of a set of integral equations. In- 
tegrability allows to introduce a one parameter deformation (the spectral parameter 
Q and study such deformed problem. One can then write down a system of func- 
tional equations, or F— system, valid for any value of ^. One can combine these 
functional equations with the knowledge of the analytic behavior of the F— functions 
in the — plane, in order to write a set of integral equations which can be solved 
iteratively, and give the desired answer. There are many directions one could try to 
follow, some of the most interesting are the following 

• It would be nice to find a physical connection between the integrable system 
that the TBA equations describe and the original integrable system. 

• It would be very interesting to extend the present construction to the full quan- 
tum problem. As a first step, one could try to compute one loop (from the 
strong coupling point of view) corrections to the above picture. This would 
allow, for instance, to distinguish between different amplitudes. 

• It would also be interesting to look for similar structures (for instance, the anal- 
ogous of the spectral parameter, etc) in perturbative computations. Related to 
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this, in |42] and subsequent papers, an operator product expansion for Wilson 
loops have been developed. This allows to use certain tools of integrability |43 1 
in order to make predictions at all values of the coupling. 

• One could hope that similar technology can be applied to related problems. 
One such problem is the computation of form factors, in which progress have 
been made recently, see [20| [44| |45|. 
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